Hole spin relaxation in p-type (111) GaAs quantum wells 
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Hole spin relaxation in p-type (111) GaAs quantum wells is investigated in the case with only 
the lowest hole subband, which is heavy-hole like in (111) GaAs/AlAs and light-hole like in (111) 
GaAs/InP quantum wells, being relevant. The subband Lowdin perturbation method is applied to 
obtain the effective Hamiltonian including the Dresselhaus and Rashba spin-orbit couplings. Under 
a proper gate voltage, the total in-plane effective magnetic field in (111) GaAs/AlAs quantum 
wells can be strongly suppressed in the whole momentum space, while the one in (111) GaAs/InP 
quantum wells can be suppressed only on a special momentum circle. The hole spin relaxation due to 
the D'yakonov-Perel' and Elliott- Yafet mechanisms is calculated by means of the fully microscopic 
kinetic spin Bloch equation approach with all the relevant scatterings explicitly included. For 
(111) GaAs/AlAs quantum wells, extremely long heavy-hole spin relaxation time (upto hundreds of 
nanoseconds) is predicted. In addition, we predict a pronounced peak in the gate- voltage dependence 
of the heavy-hole spin relaxation time due to the D'yakonov-Perel' mechanism. This peak origins 
from the suppression of the unique inhomogeneous broadening in (111) GaAs/AlAs quantum wells. 
Moreover, the Elliott- Yafet mechanism influences the spin relaxation only around the peak area due 
to the small spin mixing between the heavy and light holes in quantum wells with small well width. 
We also show the anisotropy of the spin relaxation. In (111) GaAs/InP quantum wells, a mild peak, 
similar to the case for electrons in (111) GaAs quantum wells, is also predicted in the gate-voltage 
dependence of the light-hole spin relaxation time. The contribution of the Elliott- Yafet mechanism 
is always negligible in this case. 

PACS numbers: 72.25.Rb, 71.70.Ej, 71.55.Eq, 73.21.Fg 
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I. INTRODUCTION 

Spin dynamics in semiconductor quantum wells (QWs) 
has been intensively investigated in recent years due to 
the fast development of semiconductor spintronics»i~— In 
III-V zinc-blende QW systems, spin relaxation is mainly 
determined by the D'yakonov-Perel' (DP) mechanism, 8 
which results from the momentum scattering together 
with the momentum-dependent effective magnetic field 
induced by the Dresselhaus^ and/or Rashba^ spin- 
orbit couplings (SOCs). So far, most studies fo- 
cus on the (001)ii-^ and (110) QWs^^ In these 
works, many efforts have been devoted to obtain a 
long spin relaxation time (SRT) ,iLi&igilft-g2i3rgL22i22i2* 
In (001) QWs, when the strengths of the Dresselhaus 
and Rashba terms are comparable, electron spin relax- 
ation along the [110] or [110] direction can be strongly 
suppressed. 11 ' 14 i 15 ] 18 ~ 22 i 24 ~ ^l For (110) symmetric QWs, 
the effective magnetic field due to the Dresselhaus term 
is oriented along the growth direction, which leads to 
an absent DP relaxation for electron spins along this 
direction^^i^ 

Recently, some attention has been devoted to (111) 
QWs where electron spin relaxation also shows rich 
properties ££r— Cartoixa et alM- pointed out that to 
the first order in the momentum, the effective magnetic 
field from both the Dresselhaus and Rashba terms J7i = 
(a D +a R )(k v ,-k x ,0) (with z\\ [111], x\\ [112] and y\\ [110]), 
can be suppressed to zero by setting the Rashba coeffi- 
cient or to cancel with the Dresselhaus coefficient «d 



and hence the DP spin relaxation becomes absent for all 
three spin components. However, this strict cancellation 
can not occur due to the existence of the cubic Dressel- 
haus term fl 3 = ^=[-k 2 fcj,,k 2 fe x , V2(3klk y - fej))]. Sun 

et a/.' 58 showed that with this cubic term included, the 
cancellation occurs only for the in-plane effective mag- 
netic field on a special momentum circle, which may lead 
to a peak in the density or temperature dependence of 
the SRT. Experimentally, Balocchi et alM measured the 
gate-voltage dependence of the SRT and a two-order of 
magnitude increase of the SRT can be observed, reaching 
values larger than 30 ns by applying an external electric 
field of 50 kV/cm along the growth direction of the QW. 
However, to the best of our knowledge, there are still no 
reports on hole spin relaxation in (111) zinc-blende QWs. 
In the present work, we investigate the hole spin re- 
laxation in p-type (111) GaAs QWs with only the lowest 
subband being relevant. In GaAs/AlAs QWs, the lowest 
hole subband is heavy-hole (HH) like. The Dresselhaus 
and Rashba SOC terms of the lowest HH subband $7™ 
can be written as 
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[p x (ok x ky k y ),p y (k x 3k x k y ), 
^(3k 2 x k y -k% (1) 

[a x (3k x k y — ky),a y (k x — 3k x k y ), 
a h z (3k 2 x k y - k 3 y )]eE z , (2) 



respectively, which are obtained by the subband Lowdin 
partitioning method i 41 i 42 The coefficients are given in 



Appendix [5] Here, E z is the electric field along the 
growth direction. It is noted that the Rashba term $1^ 
has not been reported in the previous literature and in- 
terestingly its component fij^ (i — x, y, z) can be strictly 
cancelled by the corresponding Dresselhaus one fij^ by 
setting the gate voltage 



E\ = -#/(«*?). 



(3) 



The cancellation gate voltages E z (about 20 kV/cm) and 
E\ (about 22kV/cm) are close to each other, indicat- 
ing that the in-plane components of the effective mag- 
netic field can be strongly suppressed. This cancellation 
leads to a strong suppression of the in-plane inhomoge- 
neous broadening 43 and hence an extremely long SRT 
for spins along the growth direction. It is further noted 
that although the out-of-plane component of the effec- 
tive magnetic field of the lowest HH subband can also 
be suppressed to zero, the cancellation gate voltage E z 
(about 10 3 kV/cm) is too large to realize. Therefore, a 
large anisotropy between the in-plane and out-of-plane 
spin relaxations is expected. 

We also investigate the hole spin relaxation in 
GaAs/InP QWs where the lowest hole subband is light- 
hole (LH) like. The Dresselhaus and Rashba SOC terms 
are then given by 



rt D = [(/^k 2 + p l2 (k 2 z ))k y , -(^k 2 + p l2 (k 2 z ))k x , 



p z (3k x ky k y )\, 



kt)}eE z , 



(4) 
(5) 



respectively. One finds that similar to the case of elec- 
trons in (111) QWSf 3 -^ the in-plane components of these 
two terms can be suppressed to zero only on a special 
momentum circle. The out-of-plane component can also 
be suppressed to zero under an unrealistic gate voltage. 
We calculate the hole spin relaxation by the kinetic 
spin Bloch equation (KSBE) approach^ 2 - with all the 
relevant scatterings explicitly included. As the Elliott- 
Yafct 44 (EY) mechanism is the leading spin relaxation 
mechanism in hole spin relaxation in bulk,— we calculate 
the SRT due to both the DP and EY mechanisms in the 
present work. For (111) GaAs/AlAs QWs, we find that 
the EY mechanism is negligible unless the DP mechanism 
can be greatly suppressed. By tuning the gate voltage, 
we predict a pronounced peak in the SRT where the SRT 
can be extremely long (upto hundreds of nanoseconds) . 
This property is favorable for spin manipulation, storage 
and device design. This peak results from the suppres- 
sion of the unique inhomogeneous broadening 43 in (111) 
QWs. The effects of hole density, temperature and impu- 
rity density on the suppression of the spin relaxation are 
investigated. We also analyze the anisotropy of the spin 
relaxation. For (111) GaAs/InP QWs, a peak is also pre- 
dicted in the gate- voltage dependence of the SRT. How- 
ever, the peak becomes less pronounced compared with 
the case in (111) GaAs/AlAs QWs. The contribution of 
the EY mechanism is demonstrated to be negligible to 
the predicted properties. 



This paper is organized as follows. In Sec. II, the ef- 
fective Hamiltonian of the lowest hole subband in (111) 
GaAs/AlAs and GaAs/InP QWs is derived. The KSBEs 
are also constructed in this section. Then in Sec. Ill, we 
investigate the hole spin relaxations in GaAs/AlAs and 
GaAs/InP QWs respectively. We summarize in Sec. IV. 



II. EFFECTIVE HAMILTONIAN AND KSBES 

We start our investigation from the p-type (111) 
GaAs/AlAs and GaAs/InP QWs. The HH and LH can 
be described by the 4x4 Hamiltonian 



H d 



Hr 
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H e + V{z)h + H E h, (6) 



since the contribution from the conduction and split- 
off bands is marginal in the present calculation. Here, 
Hl = H L + H^ is the Luttinger Hamiltonian 46 with 

H L corresponding to the part with k x ,y = 0, H$ v8v 
is the Dresselhaus SOC of the HH and LH bands^ and 
H e stands for the strain Hamiltonian based on the the- 
ory of Bir and Pikusi 42 ' 47 Their expressions are given in 
Appendix [B] V(z) is the confinement applied along the 
growth direction within the infinite-depth well potential 
approximation and I4 is a 4x4 unit matrix. He = —eE z z 
represents the electric field term, which is the source of 
the Rashba SOC^ 

By solving the Schrodinger equation of Hq = H L ' + 
V{z)Ii + H e , the subband energies for the HH and LH 
are given by 



E 



h,l 



n z 2 n 2 h 2 
2m z a A 



te'-e 



s> 



(7) 



where n z stands for the subband index along the growth 
direction, a is the well width and m'l' 1 = ™°_ rep- 



7iT 2 73 

resents the effective masses of the HH and LH along the 
z-axis. 71^3 and mo stand for the Kohn-Luttinger param- 
eters and the free electron mass, respectively. The con- 
tributions of the strain Eg and E s are given in Eqs. (|B4[) 
and (|B5[) in Appendix [B] It is noted that Eg is just an 
energy shift for all the HH and LH subbands and can 
be neglected in the present calculation. For GaAs/AlAs 
QWs, the contribution of the strain is neglected due to 
the extremely small lattice mismatch and hence the low- 
est hole subband is HH like. However, for GaAs/InP 
QWs, due to the existence of strain, the lowest hole sub- 
band is LH like. In our calculation, we focus on the case 
with only the lowest hole subband being relevant while 
other subbands have much higher energies. 

We first construct the effective Hamiltonian of the low- 
est hole subband in (111) GaAs QWs by the subband 
Lowdin partitioning metho d 41 ' 42 upto the third order of 
the in-plane momentum. The effective Hamiltonian of 
the lowest hole subband in (111) GaAs/AlAs QWs can 
be written as 



rrh 

^cff 



si + (nfe + n£) 



(8) 



Here, e£ = tt-tt, with the in-plane effective mass m^ 
given in Appendix [AJ and <x are the Pauli matrices. The 
Dresselhaus and Rashba SOC terms fi{v R \ are given in 
Eqs. (JTJ) and @. To facilitate the understanding of the ef- 
fective magnetic field from these SOC terms, we plot the 
in-plane components around an arbitrary equal-energy 
surface under the cancellation gate voltage E\ [Eq. (jHJ] in 
Fig.[U It is shown that the Dresselhaus term almost com- 
pletely cancels with the Rashba one, leading to a strong 
suppression of the total in-plane effective magnetic field. 
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FIG. 1: (Color online) Schematic of the in-plane effective 
magnetic fields of the lowest HH subband in (111) GaAs/AlAs 
QWs around an arbitrary equal-energy surface. Red (blue) 
arrows stand for in-plane components of the Dresselhaus 
(Rashba) term OJC r n and green arrows represent the sum of 
the Dresselhaus and Rashba terms fi/,, with its amplitude 
enlarged by 10 for clarity at the cancellation gate voltage E v z . 

Similarly, the effective Hamiltonian of the lowest hole 
subband in (111) GaAs/InP QWs reads 



u l — J 



(Oj, 



nfcj-f, 



in which el = ^-=r- 

k 2m; 



and m\ stands for the in-plane ef- 



fective mass with its form given in Appendix [A] The 
Dresselhaus and Rashba SOC terms ^^m) of the LH 
are given in Eqs. Q and ((3J. 

We then construct the microscopic KSBE s 5 ' 12 to study 
the hole spin relaxations in GaAs/AlAs and GaAs/InP 
QWs, respectively. The KSBEs read 5 .*! 2 . 



p\c,ircr' Pk.crcr' |coh ~r Pk,(T(r' |scat ) 



(10) 



where pk t aa' represent the density matrix elements of 
holes with the diagonal ones pk,aa = fkc {a — ±1/2) de- 
scribing the distribution functions and the off-diagonal 
ones p k ,i/2 -1/2 = Pk,-i/2 1/2 — Pk denoting the spin 
coherence. Pk, era-' |coh are the coherent terms describing 
the coherent spin precession due to the effective magnetic 
field from both the Dresselhaus and Rashba terms while 
the Hartree-Fock Coulomb interaction is neglected due to 



the small spin polarization. 5 Their expressions are given 
by 



dfk,a 



dt 

dpk 
dt 



coh 



coli 



-2a[O x (k)Imp k + n y (k)Rep k ], (11) 



_[^(k) + n !/ (k)](/ kl/2 -/ k _ 1/2 ) 

-m 2 (k) Pk . (12) 



Here, fli (i = x,y,z) are the three components of the 
effective magnetic field from both the Dresselhaus and 
Rashba terms, pk.o-cr' |scat are the scattering terms includ- 
ing the spin conserving ones, i.e., the hole-hole (V^), hole- 
impurity (U 2 ), hole-acoustic-phonon (IMq^cI 2 ), polar 
hole-longitudinal-optical-phonon (I-Mq ^! 2 ) and non- 
polar hole-optical-phonon (|Mq™ p | 2 ) scatterings with 
A = LO/TO denoting longitudinal-optical/transverse- 
optical phonon and the spin-flip ones due to the EY 
mechanism4i Their detailed expressions can be found 
in Ref. |48| (Note the scatterings in Ref. 48 are given 
in bulk and one needs to transform them to two- 
dimensional cases). The scattering matrix elements VI 2 = 



and U 2 



V fZ iUQ /e(q)] 2 |J(i<fc 



[ ~. e(q) 

with Zi being the charge number of the impurity (as- 
sumed to be 1 in our calculation)^ e(q) = 1 — 

So v q\^(^Qz)\ 2 Sk a e "'Ie stands for the screening 
under random phase approximation^ the bare Coulomb 
potential vq — =g§- with Q 2 = q 2 + q 2 ; and the form 
factor \I(iq z )\ 2 = 7r 4 sin 2 y/[y 2 (y 2 - 7r 2 ) 2 ] with y = q z a/2. 



polar 1 2 _ 
Q.LOl 



[2ne 2 w^o/(<\ 2 + ql)]{^-^ l )\I(iq z )\ 2 wit\i 



\M, 

wlo denoting the LO phonon energy and kq (Koo) be- 

j-nonp 



ing the static (optical) dielectric constant.— |Af£ on|, ' J 



3hdg 



2duj\ai 



/n\ scatterings and \Mq 



\I(iq z )\ 2 for the nonpolar HH-optical-phonon 

\I{iq z )\ 2 for the 



rnonpi2 
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2dLU A a GaAs 

nonpolar LH-optical-phonon ones with do, d, acaAs, and 
wto representing the optical deformation potential, 51 ' 
the mass density^ the lattice constant of GaAs^ and 
the TO phonon energy^ respectively. r\ = 1(5) for 
A = LO (TO). The matrix element for the hole-acoustic 
phonon scattering due to the deformation potential reads 

h ~-Q\I(iq z )\ 2 with S denoting the defor- 



Iji.rdcf 12 
|JU Q,ACl 



2dv B , 



mation potential and v s \ being the velocity of the lon- 
gitudinal sound wavei 2 ^ For the hole-acoustic phonon 
scattering due to the piezoelectric coupling, |Mq ac | 2 = 
32n 2 he 2 el 4 (3q^f [j(j fe )|2 for the longitudinal phonon 

IMP' 12 327T 2 he 2 e 2 lA 1 r/ 2 ,2 , 2 , 2 

\M. n *nl -7 ->- w OB [Qx % + % Qz + 



and 



dv sl Q 
£ Q,AC 



dv s 



Qz Qx 



HlWyl'zY 



]\I{iq z )\ 2 for the transverse phonon, 
where q' x = ^q x -^q y + ^q z , q' y = ^q x + ^ qy + ^ qz 

and q' z = — y=qx H — t=9z- fst represents the velocity of 
the transverse sound wave and ei4 is the piezoelectric 
-23. The spin mixing A kk / for the lowest HH and 



constant. 



TABLE I: Parameters used in the calculation (from Ref. |42 
unless otherwise specified) 



7i 
641 
652 


6.85 

-81.93 eVA 3 

-0.98 eVA 3 


72 


2.10 

1.47 eVA 3 

4.67 eV 


73 
651 

OGaAa 


2.90 

0.49 eVA 3 

5.65 A 


CZAlAs 


5.66 A 


ainP 


5.87 A 


Cu 


11.81 


C12 


5.32 


C44 


5.94 


d 


5.31 g/cm 3 


do 
ei4 


48eV a 

1.41 xlO 9 V/m c 

5.29xl0 3 m/s c 


ULO 


35.4 meV b 
10.8 C 
2.48xl0 3 m/s c 


OJTO 


33.2 meV b 
8.5 eV c 
12. 9 C 



a Reference _50 
''Reference 16 
c Reference 23 



LH subbands in the spin-flip scattering terms are detailed 
in Appendix ICl 



III. NUMERICAL RESULTS 

We present our results obtained by numerically solving 
the KSBEs. All the parameters used in our calculation 
are listed in Table I. The initial spin polarization is along 
the QW growth direction except otherwise specified. We 
set the initial spin polarization as 2.5 %. Note that the 
SRT is insensitive to the small spin polarization.— It is 
noted that the system is always in the strong scattering 
limit in the present investigation. 



10- 



10 1 



to- 



ur 



T=30K 

N h =3xlO U cm" 2 
N,=0 




§ 
P 



d 
v 



5 10 15 20 25 30 35 40 
E 2 (kV/cm) 

FIG. 2: (Color online) HH SRT along the z-axis t Bz and the 
inverse of the inhomogeneous broadening (Sly) of the low- 
est hole subband in (111) GaAs/AlAs QWs against the gate 
voltage E z . Red solid (blue dashed) curve represents the SRT 
due to the DP (both DP and EY) mechanism. The tempera- 
ture T = 30 K, the hole density N h = 3 x 10 11 cm~ 2 and the 
impurity density JV< = 0. Note the scale of (fly) is on the 
right-hand side of the frame. 



A. HH spin relaxation in GaAs/AlAs QWs 

Before a full investigation by numerically solving the 
KSBEs, we start from an analytical analysis of the DP 
spin relaxation of the lowest hole subband, which is HH 
like, in (111) GaAs/AlAs QWs. In the strong scattering 
limit, the evolution of spin polarization along the z-axis 
can be approximated by 51 



s, = -r* p [s z (nl + si 2 y ) - s x (si x si z 



s y (si y si z )}, (13) 



with t* standing for the effective momentum scatter- 
ing time contributed by the hole-hole, hole-impurity and 
hole-phonon scatterings. 16 i 43 The evolutions of the in- 
plane components S x and S y can be obtained by index 
permutation. It is noted that unlike (001) III-V zinc- 
blende QWs where Sl z is absent and the inhomogeneous 
broadening of the spin polarization along the z-axis is de- 
termined by (fl'l) = (Sl 2 x + Sl 2 y ), 5 for (111) QWs, Sl z and 
(Sl x Sl z ) are nonzero under the effective magnetic field 
from both the Dresselhaus and Rashba terms given in 
Eqs. (JTJ) and ([2]). This makes the evolutions more com- 
plex compared with (001) QWs, i.e., 



s z = 


-r;[s z (n 2 x + nl)~s x (a x n z )} 1 


(14) 


"I = 


-t;[s x (si 2 v + si 2 z )-s z (si x si z )}, 


(15) 


Sy = 


-t;[s v {qI + q 2 z )}. 


(16) 



The exact solutions are given in Appendix [Dj Specially 
by setting the initial spin polarization along the z-axis 
as 5 Z (0), we obtain S z = S' 2 (0)e _ */' r ^ approximately by 
considering that (SI 2 .), (SI 2 ) and (Sl x Sl z ) are much smaller 
than (Cl 2 ) in our investigation. The SRT along the z-axis 
reads 



.. = r;- l (sii)-\ 



(17) 



with (Slz) being the inhomogeneous broadening. This 
inhomogeneous broadening is unique in (111) QWs and 
has not been reported in the literature. 

We then numerically calculate the spin relaxation of 
the lowest hole subband, which is HH like, in (111) 
GaAs/AlAs QWs where the well width is chosen to be 
8 nm. We first investigate the SRT due to the DP mech- 
anism and plot the gate-voltage dependence of the SRT 
along the z-axis at T = 30 K with the hole density 
Nh = 3 x 10 11 cm~ 2 by red solid curve in Fig. [2j It 
is seen that the SRT can be effectively modulated by the 
gate voltage and a peak appears at E z m 22 kV/cm. Un- 
der this gate voltage, the SRT is as long as 283 ns which 
is about four orders of magnitude larger than the one 
without any bias. This peak results from the suppres- 
sion of the inhomogeneous broadening (SI 2 ) as shown by 
green dotted curve in Fig. [5J This suppression originates 
from the cancellation of the y-component of the effective 
magnetic field Sl y [Eq. ©]. 

Since the suppression of the DP spin relaxation is de- 
termined by the cancellation of the y-component of the 
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FIG. 3: (Color online) HH SRT along the z-axis t Sz of the 
lowest hole subband in (111) GaAs/AlAs QWs as a function of 
(a) the hole density with temperature T — 30 K and impurity 
density Ni = 0; (b) temperature with the hole density Nh = 
4 x 10 11 cm -2 and Ni = 0; and (c) impurity density with 
T = 30 K and N h - 3 x 10 11 cm' 2 . Red curves with x stand 
for the SRTs due to the DP mechanism under zero bias. Green 
curves with □ (blue ones with •) represent the SRTs due to 
the DP (both DP and EY) mechanism at the cancellation 
gate voltage E z = 22 kV/cm. (d) HH SRT along the x (y)- 
axis t 3x (T 3y ) against the gate voltage E z with T = 30 K, 
N h = 3 x 10 11 cm" 2 and Ni = 0. 



effective magnetic field Q y [Eq. Q], it is robust against 
the hole density, temperature and impurity density. In 
Fig- EJa)-(c), we plot the hole density, temperature and 
impurity density dependences of the SRT under zero bias 
and the cancellation gate voltage E z = 22 kV/cm as red 
curves with x and green ones with □, respectively. It 
is shown that the SRT at the cancellation gate voltage 
is always four orders of magnitude larger than the one 
under zero bias. In addition, we also observe peaks in 
the hole density and temperature dependences with the 
peak locations fixed under different gate voltages. The 
peak in the density dependence of the SRT has been 
theoretically predicted by Jiang and WvM- and experi- 
mentally confirmed^— recently, which is attributed to 
the crossover from the nondegenerate to degenerate limit. 
The peak in the temperature dependence of the SRT was 
predicted by Zhou et al£2- and experimentally realized by 
Leyland et al.^ Ruan et al.^ and Han et al., 5A which 
is solely caused by the Coulomb scattering^ The loca- 
tion of the peak is insensitive to the gate voltage just the 
same as the case for electrons in (111) QWsi2& 

It is further noted that the EY mechanism is usually 
the major mechanism in hole spin relaxation in bulk 
GaAs4^ In order to rule out the influence of the EY 
mechanism to the predicted phenomena, we plot in Fig. [5] 
the gate- voltage dependence of the SRT due to both the 
DP and EY mechanisms as the blue dashed curve. It 



is seen that the EY mechanism can be important only 
around the cancellation gate voltage where the DP spin 
relaxation is greatly suppressed. Away from the cancella- 
tion gate voltage, the contribution of the EY mechanism 
is negligible since the spin mixing between the HH and 
LH is extremely small due to the large energy splitting 
between the HH and LH subbands in QWs with small 
well width. Moreover, since the EY mechanism is im- 
portant at the cancellation gate voltage, we also plot in 
Fig. EJa)(c) the hole density, temperature and impurity 
density dependences of the SRT due to both the EY and 
DP mechanisms under this gate voltage (blue curves with 
• ). We find that the peaks in the hole density and tem- 
perature dependences disappear since the SRT caused 
by the EY mechanism decreases with the hole density 
and temperature^^^ Moreover, the SRT due to the EY 
mechanism also decreases with the impurity density— ^» 
However, away from the cancellation field, the spin relax- 
ation is determined by the DP mechanism only and the 
peaks still exist in the corresponding temperature and 
density dependences. 

Finally, we address the anisotropy of the spin relax- 
ation with respect to the spin polarization direction. The 
SRT along the z-axis t Sz is given in Eq. (1171) . Similarly, 
one obtains the SRTs along the x- and y-axes as 






n 2 



■W 



(18) 
(19) 



by setting the initial spin polarizations along the x- and 
y-axes correspondingly. By comparing Eq. (|17p with 
Eqs. (|T5)) and (fTi?|) , one finds a strong anisotropy between 
the out-of- and in-plane spin relaxation, since Q z is about 
one order of magnitude larger than the in-plane compo- 
nents Sl x and n y in (111) QWs with small well width. In 
addition, a marginal anisotropy also arises between the 
in-plane spin relaxations from Eqs. ([TBI and (IT51) . The 
numerical result of the SRT along the x (y)-axis r Sx (r s ) 
is shown in Fig. EJd) by red curve with x (blue one with 
□). One finds from the figure that both t Sx and t s are 
about two orders of magnitude smaller than t Sz . 



B. LH spin relaxation in GaAs/InP QWs 

We also investigate the spin relaxation of the lowest 
hole subband, which is LH like, in (111) GaAs/InP QWs 
with the well width a = 20 ran. The gate- voltage depen- 
dence of the SRT along the z-axis due to the DP mecha- 
nism under different temperatures is shown in Fig. 0] and 
a peak is also observed. We find that the peak becomes 
less pronounced than that in GaAs/AlAs QWs since the 
cancellation only occurs on a special momentum circle 
in GaAs/InP QWs. 35 In addition, when the temperature 
increases, the peak location changes, similar to the case 
for electrons in (111) GaAs QWs^ It is noted that the 
contribution of the EY mechanism is in the order of 10 ns 
and is therefore always negligible in this case. 




-160 -110 -60 -10 40 

E z (kV/cm) 

FIG. 4: (Color online) LH SRT along the z-axis t s , of the 
lowest hole subband in (111) GaAs/InP QWs against the gate 
voltage E z . Red curve with x (blue one with □) stands for 
the case with temperature T = 40 (60) K. The hole density 
Nh = 1 x 10 12 cm -2 and impurity density Ni = 0. 



IV. SUMMARY 

In summary, we have investigated the hole spin relax- 
ations in p-type (111) GaAs/AlAs and GaAs/InP QWs 
with only the lowest hole subband being relevant. The 
lowest hole subband is HH like in the former case while 
LH like in the latter one. In both systems, we uti- 
lize the subband Lowdin perturbation method to ob- 
tain the effective Hamiltonian including the Dresselhaus 
and Rashba SOCs. Under a proper gate voltage, the to- 
tal in-plane effective magnetic field in (111) GaAs/AlAs 
QWs can be strongly suppressed in the whole momen- 
tum space, while the one in (111) GaAs/InP QWs can 
only be suppressed to zero on a special momentum circle. 
Therefore, tuning the gate voltage is an efficient way to 
modulate the hole spin relaxation in (111) QWs. 

We calculate the hole spin relaxation due to the DP 
and EY mechanisms by numerically solving the KS- 
BEs with all the relevant scatterings explicitly included. 
For (111) GaAs/AlAs QWs, the contribution of the 
EY mechanism is negligible unless the DP mechanism 
is strongly suppressed. In addition, we predict a pro- 
nounced peak in the gate voltage-dependence of the HH 
SRT where an extremely long SRT (upto hundreds of 
nanoseconds) can be reached. This peak comes from 
the suppression of the unique inhomogeneous broadening 
in (111) GaAs/AlAs QWs, specifically, the cancellation 
of the y-component of the effective magnetic field. The 
suppression of the spin relaxation is robust against the 
hole density, temperature and impurity density. A strong 
anisotropy between the out-of- and in-plane spin relax- 
ations is also addressed. In (111) GaAs/InP QWs, a peak 
is also predicted in the gate- voltage dependence of the LH 
SRT. However, the peak becomes less pronounced com- 
pared with the case in (111) GaAs/AlAs QWs since the 
cancellation only occurs on a special momentum circle 
in GaAs/InP QWs. The contribution of the EY mecha- 
nism is always negligible in this case. This investigation 



suggests that the HH in (111) GaAs/AlAs QWs offers 
unique properties for the spintronic devices and calls for 
experimental investigations. 
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Appendix A: EFFECTIVE MASS AND 
SPIN-ORBIT COUPLING COEFFICIENTS 

The in-plane effective mass m^ in Eq. © and the Dres- 
selhaus (Rashba) strengths f% (aj 1 ) (i = x,y,z) of the 
lowest hole subband, which is HH like, in Eqs. ([T]) and 
(|2|) are given by 
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The in-plane effective mass m\ in Eq. ©, the Dressel- 
haus strengths l x 1 ' 12 , (5 l z and the Rashba strengths a l x z 
of the lowest hole subband, which is LH like, in Eqs. (|4| 
and (0 read 
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In above equations, 7^ are Kohn-Luttinger parameters; 
641,42,51,52 stand for the coefficients of the Dresselhaus 
SOC of the valence band; A^„ 2 (A" irl2 ) is the energy 
difference between m-th and 712-th subbands of HH (LH); 
Ajj* n is the energy difference between ni-th HH and n 2 - 
th LH subbands. B t = ^ 2 ( 2 72+73) , B 2 = ^(-^+73) 
and Bs = v^±gia) . Cl = ( 6 41 + |fe 42 _ i 652 )^, 
C 2 = ^(641+f 642-651+652), C 3 = ^^(2642+652), and 



C 4 = |(6 42 + 651 - 652). P nin2 = 4a " 1 "4-jy" 1] ( 1 - 
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Appendix B: LUTTINGER HAMILTONIAN FOR 
(lll)-ORIENTED III-V ZINC-BLENDE CRYSTAL 

The Luttinger Hamiltonian Hl, the Dresselhaus SOC 
Hg vSv and the strain Hamiltonian H t in Eq. ([5]) are given 
in the following. The Luttinger Hamiltonian Hl can be 
written as 4 ^ 
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2^-diag(7i -273, 71 +27 3 , 71 + 27 3 , 71 - 27 3 ) corresponds 
to the part of the Luttinger Hamiltonian Hl by setting 

The Dresselhaus SOC H\ v8v is given by 42 
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(i = x, y, z) represent the spin-3/2 angular momentum 
matrices and cp denotes the cyclic permutation of the 
preceding terms. 

The Bir-Pikus strain Hamiltonian H e read s 42 ' 47 
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Here, D' u and Dd stand for the deformation potential 
constants; ej, and e a denote the diagonal and off-diagonal 
components of the strain tensor; a e and a s represent the 
lattice constants of epilayer (GaAs) and substrate mate- 
rials (AlAs/InP); and C\\, C\i and C44 are the stiffness 
constants. 42,47 



Appendix C: SPIN MIXING A kk , FOR THE 
LOWEST HH AND LH SUBBANDS 



The spin mixing Ak,k' in the spin-flip scattering 
induced by the EY mechanism reads A k k / = I — 

a 2 x 2 unit matrix. For the lowest hole subband, which 
is HH like, in GaAs/AlAs QWs, S^ } is given by 
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(l|Jify|n z ) in which 1 (n z ) stands for the first 
(n z -th) subband and Hij represents the i-th row and j- 
th column of the Hamiltonian matrix i?4 X 4 in Eq. ([5]). 
A^„ is the energy difference between the first and n z -th 
HH subbands and A^ is the energy difference between 
the first HH and n z -th LH subbands. 

For the lowest hole subband, which is LH like, in 



GaAs/InP QWs, s£ ] is given by 
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where A" n is the energy difference between the first and 
n z -th LH subbands. It is noted that S k , can be obtained 
by replacing k in S* k with k'. 



Appendix D: EXACT SOLUTIONS OF Eqs. flTU) 

to dm) 



M 2w 

- O 2 ) + u)e-^ + ((tf z tf x ) + u)e-^% 

(Dl) 

5, = ?M m i-nt)+u)e-^;t + m t-tf z ) + u) 

Zu 

u 

(D2) 



S y - ^(Oje-^+^W*, 



(D3) 

With A lj2 = i((fi 2 + 2fi2 +fi 2^ ±M ) andu = (Q2^ + ^2^ 
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